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Abstract 

The four-point Green's function of static QED, where a fermion and an antifermion are 
located at fixed space positions, is calculated in covariant gauges. The bound state spec- 
trum does not display any abnormal state corresponding to excitations of the relative 
time. The equation that was established by Mugibayashi in this model and which has 
abnormal solutions does not coincide with the Bethe— Salpeter equation. Gauge transfor- 
mation from the Coulomb gauge also confirms the absence of abnormal solutions in the 
Bethe— Salpeter equation. 
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One of the puzzling features of the Bethe— Salpeter equation [I], 0, ^| is the appearance 
in the bound state spectrum of abnormal solutions corresponding to excitations of the 
relative energy (or relative time) |3|, || . These states do not have nonrelativistic counter- 
parts, neither do they appear in the hamiltonian formalism. The first example of such 
states was provided by the Wick— Cutkosky model f| |6|, where, in the one scalar photon 
exchange approximation, an infinite number of new states, with a new quantum number, 
appeared in the bound state spectrum for values of the coupling constant greater than 
a critical value. It was also established that these states did couple to the two-particle 
Green's function @. 

The fact that the presence of abnormal states is due to the excitations of the relative 
time was checked by the analysis of the problem in the static model ||, in which two 
fermions are held fixed at definite positions in space. In this case the whole set of excited 
states should come entirely from the relative time excitations. It was again found that 
in the one-photon exchange approximation an infinite number of abnormal states exist 
when the coupling constant exceeds a critical value. 

It was suspected that the appearance of abnormal solutions might be due to the 
one-photon exchange approximation and that multiphoton exchange contributions might 
remove them from the spectrum. In this connection the exact static model was analyzed 
by Mugibayashi ||. Using hamiltonian formalism he derived a second order differential 
equation that he identified with the Bethe— Salpeter equation. This equation had again 
an infinite number of abnormal states (for massless photons); the situation was even worse 
than before, since these were now appearing for any value of the coupling constant. 

The trouble with Mugibayashi's result is that it is in obvious contradiction with exper- 
imental data: no such states are observed in QED. [Mugibayashi's scalar interaction can 
be transposed to QED.] On the other hand, the one-photon approximation results might 
still escape experimental tests, since the QED coupling constant a is far below from the 
critical value at which abnormal states appear. 

These troublesome results and conclusions have led us to reanalyse the static model. 
The key ingredient in our analysis is the fact that the Green's function of the static 
model is exactly calculable and hence an explicit check of its poles provides indication 
about the bound state spectrum. Only one pole, that of the (normal) ground state, 
is found. The contradiction that exists with Mugibayashi's result stems from the fact 
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that his equation is not the Bethe— Salpeter equation, but rather a secondary equation 
of the theory. Calculating the Bethe— Salpeter kernel from Feynman diagrams or from 
the inverse of the Green's function, it is found, contrary to Mugibayashi's observation, 
that the kernel contains, in covariant gauges, an infinite series of multiphoton exchange 
diagrams. 

The above results are also verified in the Coulomb gauge, where the Bethe-Salpeter 
equation has only one bound state. Gauge invariance of the theory can then be used to 
provide additional justification to the absence of abnormal solutions of the Bethe-Salpeter 
equation in other gauges. 

We therefore conclude that the Bethe-Salpeter equation does not have any abnormal 
solution in the exact static model. Since it is unlikely that excitations of the space 
coordinates induce by themeselves relative time excitations, one should expect that this 
result remains also true in the more general four- dimensional theory. Our result confirms 
Wick's conjecture || that abnormal states are spurious solutions due to the one-photon 
exchange approximation. 

The absence of abnormal solutions in the static model has also its relevance for Heavy 
Quark Effective Theory, which is formulated as an expansion around the static model 

We now turn to the details of the calculations. We consider spinor electrodynamics 
in the static limit where one fermion and one antifermion are held fixed at definite space 
positions. The corresponding lagrangian density, in an arbitrary linear covariant gauge 
characterized by a parameter £, is: 

£ = V>i(i7o<9° - mi - gxjoA )^ + 4> 2 (i-f d° - m 2 - g^JoA )^ 

--F UU F" V - , 1 A daA^) 2 . (1) 
4 2(1-0 

Since the spacelike components of the 7-matrices are absent, the fermion fields can be 
classified according to the eigenvalues of the matrices 70 (+1 for the fermion, —1 for the 
antifermion). In the Feynman gauge, the above lagrangian density becomes equivalent, 
for its mutual interaction part, to that of the scalar interaction considered in Ref. [f|. 

The four-point Green's function describing the scattering of the two particles is defined, 
up to a normalization factor, by the functional integral: 
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G(x 1 ,x 2 ,x' 1 ,x' 2 ) = J D^D^DA^x^^^M^i^ 

= < V'i(^i)V'2(4)^ 2 (^)? 1 (xi)e^^ +i / JoAO > , (2) 

where S = J d 4 xC(x) is the action and J is an external source for the field A . To calculate 
this integral we eliminate the interaction term in S by the following transformations of 
the fermion fields: 



il>i(x) -> e- i9i ^ A ° {x) ^{x) , ^(x) -> e^^ViW (i = 1, 2) , 



(3) 



where designates a primitive of with respect to x , being defined in 

momentum space as the principal value of i/k°. The Jacobian of these transformations is 
unity and we obtain for G: 

G = < ^ 1 {x 1 )e- i9l ^ A ° {xi H 2 {x' 2 )e- i92 ^ A ° {x '^ 

Xe iS (ip,^,A)+i J Jo A > ^ 

where £0 is the free action . 

The remaining integration in the functional integral is gaussian and yields in a straight- 
forward way: 



where G is the free two-particle propagator, 

Co = G10G20 1 

G i0 ( Xi - x' t ) = 9(x° - x^e-^-^S 3 ^ - xj) (i = 1,2), 



(5) 



and J° is: 



J (/c, x±, x 2 , X-^, X2) — i~. 

Kq 



g!(e ik - Xl - e ik -<) - g 2 (e ih - X2 - e ifc ^) 



(6) 



(7) 



For J° = one has: 

G = G e <^ 



91 (e ikl l - e ik - x l)-g 2 (e ik - x 2 - e ik ' X 2 ) 



(fc 2 +«)fcg 



(8) 



In the argument of this exponential the terms proportional to gf (i = 1, 2) contribute to the 
mass and wave function renormalizations. After carrying out the latter renormalizations 
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(we shall continue denoting by the same notations the renormalized masses) and setting 



9i = 92 = g 



a 



9_ 

Air 



4vr 2 



(9) 



we obtain from Eq. (H) for the renormalized Green's function: 



G = G e 



Xa(xi-x 2 )-xt(xi-x' 2 )-xz(x2-x' 1 )+X(;(.x' 1 -x' 2 ) 



xe 



h^ixO-x'^+^^xO-x'V) 



(10) 



where \£, an d h,^ (i = 1, 2) are defined as: 



d A k 



,ik.x 



(2tt) 4 k 2 + ie 



x° + it 



7-2 
A 



-ikox 



1) 



-iX 



2r 
(C + 2) 



/c 2 + ze (-A; + ie)(k + ze) 
+ (i+2) 1 / z£ !+i£^ 



r — x u 



ln(— mf(x 



02 



rg)+ze) 



1,2) 



(11) 
(12) 



where r (— 0) is the renormalization short distance regulator. Notice that x is an even 
function of x°. The functions come from the one-fermion propagator renormalizations. 
They represent the cloud of scalar and longitudinal photons usually appearing in covari- 
ant gauges ||11|| . Their effect in momentum space consists in replacing the pole of the 
propagator by a power law singularity. A similar effect is also contained in the second 
term of • 

In the absence of radiative corrections (h^ = 0), formula ( ]T0"| ) can also be obtained by 
directly summing the Feynman diagrams corresponding to the multiphoton excahanges 
and contributing to the fermion-antifermion scattering amplitude. Because of the linear- 
ity of the denominators of the fermion propagators in momentum space [Eqs. @], the 
summation technique of the eikonal approximation [|l^] can be applied in exact fashion, 
with the fermions kept off their mass shell. When radiative corrections are present, the 
calculations are more involved. The emission amplitude of n photons from a fermion 
line is completely determined, through successive uses of Ward-Takahashi identities, in 
terms of the fermion propagator. After this step, the eikonal summation technique can 
be applied with appropriate adaptations and formula ( ]TD| ) is reached. 
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In the following we shall use for the total and relative variables the notations: 

^ _ X\ + X 2 _ x[ + x 2 ^ _ _ j£, 

P = Pi+p 2 , P = ^(pi-P2) ■ (13) 

To extract from the Green's function (|H]) the bound state spectrum of the theory we let 
the time t tend to infinity, keeping x° and x'° fixed. Because of the cluster decomposition 
H, the Green's function behaves in general for infinite time separation as: 



G(x 1 ,x 2 ,x[,x' 2 ) ~ ^0 n (x 1 ,x 2 )0 n (x' 1 ,x / 2 ) = Y,M x )e iP ° nt <t> n { x ') ■ ( 14 ) 
In this limit the Green's function (|T0|) yields only one exponential factor and behaves as: 

G ~ e -i(mx+m2-a/r)t _ ^ 
t^oo 

This behavior signals the presence of one bound state with energy 

a 

P = mi + m 2 , (16) 

r 

which is gauge invariant. Notice that the contributions of the scalar and longitudinal 
photon clouds in the functions \i an d have cancelled each other in the limit t — > oo. 
This is a consequence of the electric neutrality of the bound state considered here. 

The above result about the bound state spectrum can also be obtained by analyzing 
the Green's function in total momentum space. After the isolation of the ground state 
pole fllBD , no other infinite type singularities appear. 

We conclude from this calculation that the static model has only one bound state, the 
normal one, and no abnormal states exist. 

We next turn to the analysis of the problem by means of the Bethe-Salpeter equation. 
By functional methods [TB[] one establishes the following equation for the Green's function: 

6_ 

'SJ (xi) 



(Hiod xw -m 1 )G(x 1 , x 2 , x[, x' 2 ) = i5 4 (x 1 -x' 1 )G 20 (x 2 -x 2 )-ig'y 10 —- — rG(x 1 ,x 2 ,x' 1 ,x 2 ) , 



(17) 

and a similar one with the interchanges 1 «-> 2. Here Gj (i = 1, 2) is the free propagator of 
fermion i [Eq.(^j)]. The complete expression of the Green's function [Eq. @] can then be 
used to transform these equations into integro-differential equations for G. Equivalently, 
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one can directly apply the operators G~q (i = 1,2) on the expression (|TDJ) of G. One 
obtains the following two differential equations (we omit radiative corrections, which 
actually disappear from the subsequently derived wave equations): 



(idxio ~ m i)G(xi, x 2 , x[, x' 2 ) = i5 4 (xi — ^ / 1 )G 2 o(^2 — x' 2 ) — 



xG(xi, x 2 , x[, x 



2) > 



X , i (x 1 -x 2 ) - xfai -4) 

(18a 



(id x2 o ~ m 2 )G(x 1 , x 2 , x[, x' 2 ) = i5 A (x 2 - x' 2 )G 10 (x 1 - x[) - 



X^X-2 ~ Xi) - x^{x 2 - x Xi 



xG(xi,x 2 , x[,x 



2) i 



;isb) 



where the prime on x designates the derivative with respect to the temporal argument: 





x 



dx° 



xd x ) ■ 



(19) 



x . 



[x^(a;) is an odd function of 

Continuing the procedure, one also obtains a second order differential equation for G: 

(id xl0 - mi)(id w20 - m 2 )G(x u x 2l x' 1 ,x 2 ) = i 2 8 A {x x - x[)5 4 (x 2 - x' 2 ) + I ix'[(xi ~ x 2 ) 



+ 



(x't(xi ~ x 2 ) - x'^xi - 4)) (x't(%2 - xi) - x'^(x 2 - x[) 



G(xi 1 x 2l x' 1) x' 2 ) . (20) 



Wave equations satisfied by the Bethe— Salpeter wave function are obtained by for- 
mally taking in the previous equations for G the limit t — > oo [Eq. (|T3|)1, using the cluster 
decomposition (|14]) and eliminating by integration on X'° all the amplitudes <p n but one 
§. One finds: 

Att" 



(i<9io - mi)0(xi,x 2 ) 
(id 20 - m 2 )(j)(x 1 ,x 2 ) 

{id w - m x ){id 2 Q - m 2 )0(xi,x 2 ) 



2r 



40) + 



Att 
2r 



[Xi,X 2 ) , 



(Xl,X 2 ) , 



(21a) 
(21b) 



ix'[( x ) ~ (xf(x) - (^) 2 



{x 1 ,x 2 ) • 



(22) 



In the Feynman gauge these equations reduce to those obtained by Mugibayashi |J with 
the hamiltonian formalism. 

Equations ( |21a| )-( pT5D yield for <p a single bound state with the energy ([16]) and wave 
function 

= Ce -i(mi-rn 2 )xV2 e -i(m 1+ rn 2 -a/r)X° ^ } ( 2 3) 



with C, a constant. 

Equation (|22~D is best analyzed in euclidean space for the relative time variable. Using 
Wick rotation || [| and then setting p = ip$ or x° = —ir, Eq. (p2|) becomes (for 
simplicity we consider the equal mass case mi = m 2 = m): 



-T~ m)0= -^- 2 (r 2 + r^ 

2 



(2 + Or 2 + (2 -Or 2 



2r/ 



iarctan(^ ir 



1 2 



• (24) 



r r 2(r 2 + r 2 ) 

For large \t\ the potential behaves as —(1 + £/2)A 2 7r/(2r|r|), indicating that for £ > —2 
there are an infinite number of normalizable solutions for any value of A. These additional 
solutions were identified with the abnormal solutions of the Bethe— Salpeter equation. On 
the other hand, had we kept in the above equation only the linear term in A (coming from 
the one-photon exchange contribution), we would have faced the same situation as in the 
Wick— Cutkosky model |5], or its static analogue of Ref. ||. In this case the additional 
solutions exist only for A > 1/(4(1 + ^/2)). 

It is clear that these abnormal solutions do not appear as poles of the Green's function 
because they do not satisfy the first order equations ( j21aj) - (|21b| ). At this stage the main 
objection one could raise is the fact that Eqs. (18a)- (|18b|) are not fundamental equations 
of Quantum Field Theory and are specific to the static model; they would not survive 
in more general cases and hence only Eq. ( [22] ) would remain as a bound state equation, 
leading to the existence of abnormal solutions. Leaving aside for the moment the feature 
that the very existence and the energies of the abnormal solutions are manifestly gauge 



dependent and in any event the latter could not appear as poles of Green's functions ||14 



Eq. fl2"2|) suffers from the following main drawback: in spite of appearances, it is not the 



Bethe— Salpeter equation. The reason for this is that the potential in Eq. (22) does not 
represent the Bethe— Salpeter kernel. If this was the case, then we would conclude that 
the series of multiphoton exchanges in the kernel stops at two photon exchanges with a 
corresponding local expression in x-space. The contribution of the two-photon exchange 
diagram to the kernel can be explicitly calculated. In the Feynman gauge it reads: 



K 2 { Xl , x 2 , x[, x' 2 ) = \ 2 6(x\ - x?)6{x 2 - x 2 °) 



-im\ (x®—x'®)—im2 {x^—x'^) 



\(x l - x' 2 ) 2 - ie][(x 2 - x[) 2 - it] 
x <J 3 (xi-xi) 5 3 (x 2 -x' 2 ) . (25) 
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This is not a local operator in the temporal variables and does not cope with the 0(A 2 ) 
terms of the potential of Eq. 



More generally, the 0(X ) terms of Eq. fl20|) , which are multiplicative functions and 



act on the four arguments of the Green's function, do not represent the explicit expression 
of the Bethe— Salpeter kernel, which should act on two arguments only. The latter can 
be calculated either from the Feynman diagrams or iteratively by identifying its global 
action with the multiplicative functions of Eq. (p0|). The general expression of the kernel 
cannot, however, be put in a compact form and will not be reported here. Its main 
feature is that the series of multiphoton exchange diagrams does not stop at a finite order. 
Because the Bethe— Salpeter kernel is directly related to the inverse of the Green's function 
[K = Gq 1 — G^ 1 ], the bound state spectrum of admissible states of the Bethe— Salpeter 
equation should be the same as that of the Green's function itself. 

The absence of abnormal states in the Bethe— Salpeter equation can also be verified 
using gauge invariance of the lagrangian density (0) (up to the gauge-fixing term) and 
working in the Coulomb gauge. In this gauge, the expression of the Green's function is 
obtained from Eqs. (|lT)|)-(|TT|) by the replacement of the photon propagator contribution 
by — 1/k 2 . This amounts to replacing by zero and %g D Y Xc with: 

Xc(x) = ~\x°\ . (26) 
2r 

For large time separations, the Green's function behaves as: 

G C ~ e -i(mx+m 2 -a/r)t . ^ 
t— >oo 

it displays one bound state with energy given by Eq. fllBT). 

The equivalent equations to Eqs. ( |21aj ), ( |21b| ) and (|22| ) become: 



(z<9io - mi)<f>c(xi, x 2 ) = -~9{-x a )<l>c{x u x 2 ) , (28a) 

r 

{idw - m 2 )(t)c{xi ) x 2 ) = O(x )(j)a{xi,x 2 ) , (28b) 

r 

(id 10 -m 1 )(id 20 -m 2 )(f)c(xi,X2) = i-8(x )4>c(xi,x 2 ) . (29) 

r 

The latter equation is actually the Bethe— Salpeter equation in the Coulomb gauge, be- 
cause the instantaneity of the interaction in this gauge reduces the Bethe— Salpeter kernel 
to the one-photon exchange diagram. This equation has clearly one bound state solution, 
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independently of Eqs. fl28ap and ( |28b| ), which confirms the fact the Bethe— Salpeter equa- 
tion should by itself contain the full information about the bound state spectrum. 

The gauge transformation operator between two covariant gauges, or a covariant gauge 
and the Coulomb gauge, for the Green's function is obtained by taking the ratio (in x- 
space) of their corresponding expressions; it reduces to the ratio of the exponential factors 
with the x an d h functions. Taking the large time separation one also deduces the gauge 
transformation between the corresponding Bethe— Salpeter wave functions. One finds: 



z i[x ^ x) - x ^ x)] ^(x 1 ,x 2 ) , <pc(x 1} x 2 ) =e i[xc{x) ^ {x)] <p ! .(x 1 ,x 2 ) . (30) 



The above transformation laws, which reflect general properties of the Bethe-Salpeter 
wave function |L5|| , imply that if a bound state exists in one gauge, it should also exist in 



any other gauge. This results from the fact that the exponential factors involve the differ- 
ence of two x's and hence do not qualitatively modify (in euclidean space) the dominant 
asymptotic behavior of the wave functions. 

It can be checked that these transformations applied on Eqs. fl21a|) - (|2TED yield the cor- 
responding equations in the new gauge. However, the equivalence of the Bethe— Salpeter 
equations of two different gauges cannot be established in strong form. A general feature 
of the Bethe— Salpeter equation is that it is only weakly invariant under gauge trans- 
formations, that is that its (infinitesimal) variation is proportional to the starting wave 
equation. [This behavior can be established by writing the equation in the form G~ l <p = 
and then using the infinitesimal gauge transformations of G and 0.] In particular, using 
in Eq. (|29| ) transformation ([30]) one obtains an equation that is weakly equivalent to the 
Bethe-Salpeter equation of the covariant gauge. Nevertheless, the general property of the 
invariance of the bound state spectra under gauge transformations |14| , explicitly verified 
on transformations fl30"|) , ensures us that the latter equation has also one bound state. 

The previous calculations can also be repeated with a massive photon with mass \x. 
For large time separations, the Green's function behaves as: 

Q ^ g-i(mi+m 2 -ae"'"'/r)( (31) 

We deduce that there is only one bound state, the normal one. 

In summary, the exact expression of the four-point Green's function in the static 
model does not display any abnormal state. In the one-photon exchange approxima- 
tion (in covariant gauges), the situation is similar to that of the Wick— Cutkosky model, 
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where abnormal solutions appear for sufficiently large values of the coupling constant. 
The present study shows that the multiphoton exchange contributions sum up to sweep 
away the previous phenomenon and to reestablish a normal structure of the bound state 
spectrum. Since the phenomenon of the abnormal states is solely due to the relative time 
excitations, the results obtained in the static model, where spacelike motion is frozen, 
should survive in the more general case of four- dimensional theory. 
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